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ABSTRACT 


Borders are monitored by a variety of moving and stationary sensors, e.g., patrol 
agents, video cameras, ground sensors, UAVs, etc. This paper introduces a model for a 
moving sensor that patrols a perimeter that is infiltrated by malevolent agents (targets). 
Targets arrive according to a Poisson process along the perimeter with a certain 
distribution of arrival location, and disappear (renege) a random amount of time after 
their arrival. The measures of effectiveness (MOEs) presented in this paper are the target 
detection rate and the time elapsed from target arrival to its detection (waiting time). We 


study two types of sensor trajectories that are periodic and maintain constant speed: 


1. The sensor moves from a starting point to a certain location and then leaps 
instantaneously back to the starting point. 


2. The sensor moves back and forth between two points. 


The controlled parameters (decision variables) are the beginning and end points of 
the patrolled sector. Properties of these trajectories are demonstrated in great generality. 
The results give decision makers a powerful tool for optimally deploying and operating a 


variety of sensors in an area of interest. 
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EXECUTIVE SUMMARY 


Borders are monitored and protected by border patrol agents, video cameras, 
ground sensors, physical barriers, land vehicles and manned aircraft. The diverse nature 
of U.S. border defense is challenged by an equally diverse array of threats, ranging from 


terrorists to drug smugglers, arms dealers and human traffickers. 


This paper presents a model for a moving sensor patrolling a perimeter. The goal 
is to find a patrol policy for the sensor that minimizes the infiltration across the perimeter 
(maximize detections). We assume the sensor follows a periodic trajectory at constant 
speed. The controlled parameters (decision variables) of the trajectory are the beginning 


and end points of the patrolled sector. 


Theoretical framework is established for two general trajectories: a Back and 


forth (B/F) trajectory and a Leap to origin (LTO) trajectory. 


The LTO trajectory represents cases where the perimeter is scanned at a constant 
speed (or constant angular speed in case of a stationary camera) from the origin to a 


certain point, and then instantaneously starts to be scanned from the origin again. 


The two main assumptions in this paper are that the target arrival process is a 
constant rate Poisson process and the target arrival locations, arrival times and reneging 
times are independent. For convenience, it is assumed that the sensor’s speed is 1. It is 


then shown that the results can be adapted for any sensor speed. 


Properties of these trajectories are demonstrated in great generality. The results 
give decision makers a powerful tool for optimally deploying and operating a variety of 


sensors in an area of interest. 
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I. INTRODUCTION 


A. GENERAL 

Borders are monitored and protected by border patrol agents, video cameras, 
ground sensors, physical barriers, land vehicles and manned aircraft. The diverse nature 
of U.S. border defense is challenged by an equally diverse array of threats, ranging from 
terrorists to drug smugglers, arms dealers and human traffickers. An increasingly 
common approach to border surveillance is to use a combination of long range sensors 
and short range sensors [1]. Generally, there is a tradeoff between the range of the sensor 


and its resolution; long range sensors provide less resolution than short range sensors [1]. 


Although the moving sensor model developed in this paper can be applied to a 
variety of monitoring systems, we focus on a UAV-mounted sensor and a scanning 


camera system. 


There are two different types of UAVs: drones and remotely piloted vehicles 
(RPVs). Both are pilot-less, but drones are programmed for autonomous flight, whereas 
RPVs are actively flown, remotely, by a ground control operator. UAVs have played key 
roles in recent conflicts [2], providing reconnaissance, surveillance, target acquisition, 
search and rescue, battle damage assessments and attack capabilities. An example of such 
a system that is increasingly present in U.S. border surveillance activities is the Predator 
B UAV system, which can provide long-endurance surveillance and communications 


relay [3]. 


While UAVs are used to extend the surveillance range, especially in the absence 
of ground infrastructure, ground-based electro-optical sensors achieve greater resolution. 
A variety of electro-optical imaging sensors are employed in border surveillance. A few 
examples of the technologies in use are high-resolution imaging, motion detection, 
temperature-differentiation and night-vision devices. In addition there is wide use of 
commercially available CCTV cameras adapted for daytime surveillance, infrared video 
(IR) detection systems and laser illumination systems that enable high-resolution 


imaging. Other optical components include computer-operated pan/tilt/zoom cameras, 


visible or near-infrared illuminators for night vision with conventional cameras and 


image-intensifiers for long-range night vision [1]. 


The combination of different sensors mitigates the shortcomings of each single 
sensor. For example, IR provides detection in atmospheric conditions where video 
imaging is ineffective, even though the resolution it provides is relatively low. Laser 
illumination sensors provide a combination of long range and high resolution images [5]. 
B. PREVIOUS WORK 

Previous work on linear patrols mainly concerns anti-submarine warfare. The 
classic WWII “Search and screening” by Kooperman [7] lays the foundations for 
analyzing the performance of barrier patrols. In previous models of barrier patrols the 
sensor’s range is taken into account. Washburn, [9] and [10], supports and enhances the 
basic barrier patrol model. More work on barrier patrol had been done specifically for 
Anti-Submarine Warfare. In this work the acoustic degradation of the sonar sensor 
caused by increasing the patroller speed is incorporated into the model, and the optimal 


patrol speed is found. 


A game analysis of choice of speeds by infiltrator and patroller originated in 
Kettelle and Wagner [11], a precursor to [8], and was developed analytically for a 
continuum of speed choices by Langford [12]. 

C. THE MODEL FRAMEWORK 

This paper sets forth a model for a moving sensor that patrols a perimeter. The 
goal is to find a patrol policy that minimizes the infiltration across the barrier (maximize 
detections) when the sensor follows a periodic trajectory at constant speed. The 
controlled parameters (decision variables) of the trajectory are the beginning and end 


points of the patrolled sector. 
We now describe the two main components of the problem: sensor and targets 
Sensors. We consider two types of sensors: 


A UAV-mounted sensor that travels along a uni-dimensional perimeter, possibly 
of infinite length, and performs a back and forth motion with a constant speed as 
described in the picture below. This trajectory will be referred to as a B/F trajectory; see 


Figure 1. 


4 Arrival 
o Reneging 
* Detection 





Figure 1. _ Sensor and targets location as a function of time in a B/F trajectory. 


1. A human operated camera that scans a sector. The camera starts at the 
origin, moves at constant speed to a certain destination point, and then leaps to the origin 


instantaneously. This is referred to as a “Leap To Origin” (LTO) trajectory; see Figure 2. 
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Figure 2. Sensor and targets location as a function of time in an LTO trajectory. 


In both B/F and LTO we assume the origin and destination are constant in time 
and the speed is also constant and equals 1. Cycle time is the time it takes the sensor to 


depart from the origin and arrive back to the origin. 


Targets: Targets arrive to the perimeter according to a Poisson process with rate 
a(t), t> 0. More specifically, let T=(T; : i= 0) be the collection of arrival epochs, then the 


number of arrivals by the time t, t > 0 is given by 
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N,()= SF, <2), (1) 


where 


EN ,(t) = {ovw)du : (2) 


The arrival location of the targets in the perimeter is described as the collection of 
IID rv’s X=(X; : 1= 0) defined on the same probability space as T, and are conditionally 


independent of arrival epochs, so that X; may be dependent on T; 
P(X, €dx|T.,...,T,) = P(X, €dx|T,). (3) 


After a certain amount of time a target may no longer be available for detection; 
this can represent the fact that the target infiltrated across the perimeter successfully. We 
call the period of time elapsed from target arrival until infiltration the reneging time. The 
reneging behavior is described by a collection of IID rv’s R=(R; : I = 0) defined on the 


same probability space as T. The Ris are also conditionally independent 
P(X, €dx|T),...,T,,X,,...,.X;) = P(X, €dx|T,,X,). (4) 


The interpretation is that reneging behavior may be dependent on the arrival 
epoch and on the arrival location. 


The main vehicle for our formulation is the process N,(¢), that counts the 
number of detections by time t. The expression for N,(t) and EN,(t) depends on the 


trajectory policy, and will be developed in Sections Two and Three. 
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Examples of target behavior are shown in Figures | and 2. The targets’ arrival 
times and arrival locations are marked with triangles and the reneging time intervals are 
marked with a broken black line. If a target reneges before the next arrival of the sensor 
to its location, it exits the system without being detected and the reneging time is marked 
with a circle. If the sensor visits the target’s location before it reneges, the target is 


detected and the detection time is marked with a star. 


This paper establishes a framework for the employment of UAV and camera- 


based sensor platforms; specific contributions include the following: 


ie In the LTO setting, the measures of effectiveness (MOEs) are the target 
detection rate and the time elapsed from target arrival to its detection (waiting time). We 
obtain the optimal trajectory length and its location under general conditions. This is the 


subject of Section One. 


2; In the B/F setting, the measure of effectiveness (MOEs) is the target 
detection rate. We obtain the optimal trajectory lengths and location, and compare the 


B/F and LTO performance. This is the subject of Section Three. 


3. We study the case where detection is not instantaneous and each target is 


inspected by the sensor for a certain amount of time. 
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Ul. THE LTO TRAJECTORY 


A. THE RATE OF DETECTION AS A FUNCTION OF SECTOR LENGTH 
Our goal is to find the scanning policy that maximizes the expected detection rate, 
prior to reneging, per-cycle. It is assumed that the sensor always travels with a constant 
speed. For simplicity, the sensor speed in the following formulations is 1. In Section 
Three, we discuss the case where the sensor has a different speed than 1. We study the 
simplest possible case, where T, X and R are independent and independent of the arrival 
process, which is a homogenous Poisson process with rate 1. Additionally, we assume 


that F(-) is continuously differentiable, and F,(-) is continuous on (0,0), where F’(-) is 


the CDF of X, and F,(-) is the CDF of R. 


This assumption guarantees that F'(-) is sufficiently smooth for our optimization 


problems, and is satisfied by any rv with a continuous density function. Let 


A(u):=sup{F.(a+u)-F.(a)} Vuru>0, (5) 


and a =a(u) be the maximizer in the last equation, i.e., the best UAV origin for a cycle 
of length u. Clearly A(-) is continuously differentiable. Observe that /A(0)=0, 
h(u) >lasu—oo and hA(u) is non-decreasing. Therefore, /(-) is a distribution function 


with domain [0, ©); in particular, /(-) has a continuous density h'(-). 


Example 1. When X ~ N(0,o7) we get 


h(u) = O(u/20)—-®(-u/2c) =1-—20(-u/20), (6) 


where () is the distribution function of a N(0,1) rv, and a =-w/2. 


For a sensor with origin @ and destination @ +4, we have 


N,((k+)u) = N,(ku)+ > Ia’ SX, <a" +u,X,+(k-NusT, < X,+ku,R, >X,+ku-T,) 
i=l 
(7) 


for «2! is the number of detections at the end of the (k+1)’th cycle, and 


N,(u)=> l(a’ < X, <a" t+u,T. < X,,R, > X,-T) (8) 


i=l 


is the number of detections at the end of the first cycle. Proposition 4.4.1 of [6] implies 


that, along a location X = x arbitrary, the number of detections in cyclek,k = 2, of length 
u is anon-homogenous Poisson process with instantaneous rate P(R>1t),0<t<u. 


Using Proposition 4.4.1 of [6] once more to account for the randomness in arrival 


locations, the expected number of detections per cycle of length u is 
E(N,((k+1)u)— N,(ku)) =h(u)| Fe(tdt, (9) 
0 
where Fr(.) =1-F 2(:), Lherefore, the expected detection rate prior to reneging is 


(Fe (t)dt 
g(u):= ie) (10) 


The fact that u = EN (u) leads to a second interpretation of Equation (5): g(u) is 
the ratio of the expected number of detections by time u and the expected number of 
arrivals by time u. Accordingly, our goal in the LTO setting is to find a cycle length wu" 
that maximizes g(.). 

Remark: Changing the arrival rate of the Poisson process changes g(.) by a 


constant factor, but does not changew . 

We define g(0)=0. To complete the notation, let 7:=inf{v:F,(v)=1}, and 
i :=inf {v: h(v)=1}; both F and w can possibly be infinite. 

The following examples illustrate the scenario; we start by analyzing the case 
when the mass of R is concentrated in a finite time interval. 


Example 2. Suppose that 0<7 <0, so that R has no mass to the right of 7. 
Then Equation (5) becomes 


[Fe@at 
OO < r 
g(u) = h(u) 5 for0<u<r (11) 
he foru>r 
u 


When is it optimal to stop scanning before the cycle length reaches 7? The 


following lemma addresses this question. 
Lemma 1. Suppose 0 <7 <oo and that uh'(u)<h(u) foru>7 .Then u’ <7. 


Proof outline. It is sufficient to show that 7 = arg max,.. h(u)/u . Observe that 


u2r 


h(u)/u>O0asu—>o, and that h(r)/r>0. Now uwh'(u)<h(u)foru=r implies 


(h(u)/u)'< 0 in (7,0), so that there is no stationary point in (7,0). 


The borderline case isuh'(u) =h(u) foru>7 , and we must have hA'(u) constant 


(i.e., uniform) on u =>7 ; see Example 3 below. 


Under the conditions of Lemma | it is never beneficial to patrol a sector longer 
than r because the detections gained by covering more area will not out-weight the extra 
missed detections due to reneging. 


Two special cases of Lemma | are presented. 


Example 2.a) R is deterministic, 1.e., P(R=7r)=1 for some constant 7 >0. 


Equation (10) becomes 


h(u) for 0<u<r 


gtu)= (12) 


h(u)— foru =r 
u 


Because g(.) is non-decreasing on0<u <7, we have uw’ =argmax,...{h(u)/u}. 


If A(.) meets the assumptions of Lemma 1, thenu” =7. 


Example 2.b) If R ~ U(0,7-) , Equation (10) becomes 


nen(1- 2) for 0<u<r 
2r 


g(u)= (13) 


A(u)—— foru=r 
2u 


The last example can also be interpreted in terms of average waiting time. More 
precisely, suppose we wish to find the optimal cycle length in order to minimize the 


targets’ average waiting time. By time ¢ >w, total waiting time is given by 


10 


2 2 


=a t t u 
expected waiting time number of cycles expected waiting 
of customers that fall by time t time per cycle 


outside the sensor's 
sweeping region 


It is easy to see that minimizing the last expression is akin to 
maximizing h(u)(t—u), which is the top branch in the RHS of Equation (11). The next 
lemma addresses this scenario, but first some notation must be introduced: Suppose that 


x(t), y(t) are real functions, we say that x(t) ~ y(t) if lim, x(t)/y(@) =1. 


Lemma 2. Assume _ that 4h(.) ~~ is _ strictly concave. Then 


u (t)=arg max, ., h(u)(t—u) is unique. Moreover, if log(h'(u)) ~ —I(u), for (.) a strictly 


increasing function such that liminf /(u)/u* >1,«>0, then u(t) ~ J" (log(t)) ; 

Proof outline. For  /(u) strictly concave on (0,t) implies A(u)(t—u) _ strictly 
concave. Additionally, (h(u)(t—u))'> 0 for all wu € (0, (4) and 
(h(u)(t—u))'<0 for allu > &(t). Because (h(w)(t—-u))' is continuous and _ strictly 
decreasing, we must have (a(u" (t)\(t-—u ()) ‘0 for some 
u (t)e(e(),E(0)), andu'(¢)unique. For the last _ part,(A(w)(t—w))'= Oiff 
log h(u) —logd —u/t) = logt+log (h (u)) . An expansion in series now yields the result. 


The last lemma says that if the density h'(.) decays exponentially fast with rate 
I(.) (which is the case for most densities of practical interest), then the optimal cycle 


length (in terms of waiting time) grows in time like J“ (log(t)). 

Example 1 (continued). In the setting of Example 1, we get h'(u) = ¢(-u/2)/o, 
so that J(u) =(u/2o) and u(t) ~ 20. Jlog(t) ; 

Our next lemma considers the case of w finite. 


Lemma 3. Suppose that 0<<0o.Then uw’ <7. 


11 


Proof outline. The fact that } Fr(t)dt/u is non-increasing in u, and that 
0 


h(u) =1 for all u>a, implies that u <7. 


Lemma 3 is the scenario corresponding to the rv X taking values in some bounded 
set. It says that it is not beneficial to patrol further than w, where there are no targets 


arriving. 


Example 3. Suppose X ~ U[0,u], then Equation (1) is 


[Fe@at 
?____, for 0<u<a 
gua), " 3 (15) 


} F a(t)dt 


, foru>u 
u 


Clearly, wu’ =u. In particular, if 7 <w, any u’ € [7,7] is optimal. 
B. OPTIMALITY CONDITIONS 


We now discuss sufficient conditions for optimality in a general LTO setting. Let 


h'(u) x Fru) 1 . (16) 





V(u) = h - 
(w) | F r(t)dt 


and observe that v(.) is continuous on (0,0). Setting g'(u)=0 in Equation (1) and 
rearranging yields v(u)=0. We are interested in finding conditions that guarantee that 
there exists a unique point wu’ such that v(u’)=0, which then would imply that w’ is the 


best cycle length (because u’ =0 is not optimal, and g(w) > 0). Clearly, v(u)>0 for all 


12 


u €(0,€),€ >0, and observe that lim sup v(w) < 0. Now the continuity of v(w) ensures that 


u is unique if v'(w) <0. We summarize this discussion in the next lemma. 


Lemma 4. If v'(u)<0 for all u>0, then there exists a unique finite u that 


maximizes g(.) (and is the root of v(u') =0). 


Lemma 4 says that we increase cycle length until the point where the marginal 
increase in target detections equals the marginal decrease in lost target detections due to 


reneging. 
Remark. Dropping the condition v'(u) <0 raises the possibility of having more 


than one stationary point, in which case we need to compare all the roots of v(.). 
As an application of the last lemma, we have the following examples. 


Example 4. Suppose X ~ N(0,1) and R ~ exp(1), then Equation (13) yields u*= 
2.0481; see Figure 3. 


The positive and negative expressions in eq. (5) as afunction of u for sigma=1 tao=1 
5 
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Figure 3. _ Derivative of the rate function when X~Normal R~Exponential 


The reneging distribution being exponentially distributed leads to the same 
solution that arises in the following setting: Upon arrival, each target emits a signal of 


constant amplitude that decays exponentially in time with rate 8; this is known as a shot 
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noise process ([6], p.326). In particular, the signal strength of arrivals detected in the 


(k+1)’th cycle is given by 


Ye Mea" < X, <a" +u,X,+(k-YusT <X,+ku). (17) 
i=l 


Following the approach taken in [6], p. 326, we find the expected average signal 
strength of targets detected in an arbitrary cycle is h(u) [, e“dt/u, which is g(u) when 


R ~ exp(@). In other words, a system where targets that arrive according to a shot noise 


process, have the same solution as the system where targets renege according to an 


exponential distribution. 


Example 5. Let X ~ N(0,1) and R~U(0,2). Equation (13) yields u*=1.8626; 


see Figure 4 below. 


The positive and negative expressions in eq. (5) as afunction of u for sigma=1 ,r=2 
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Figure 4. _ Derivative of the rate function when X~Normal R~Uniform 
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Cc; CHANGES IN VELOCITY 


A change in the sensor’s velocity from 1 to v>0 changes Equation 1 to 


uly 


g,(u) = h(u) | F..(t)dt (u/v). Changing variables, the latter equals h,(u') } F.(tdt/u', 
0 0 


where h,(u'):= h(u'v). Because for v > 0 arbitrary, h,(u')is a distribution, the analysis of 
this section applies to the general velocity setting, with /,(.) in place of h(.). It can easily 
be seen that uw, <u’ and g,(u,) < g(w ) if and only if v <1; in other words, decreasing the 


velocity reduces the optimal cycle length and the optimal detection rate. 
D. SCANNING OF DISJOINT SECTORS 

Our definition of h(.) (cf. Equation (1)) constrains the camera operator to scan a 
continuous sector. In reality, however, there may exist sectors that are impassable, in 
which case the camera operator skips the forbidden sector to scan a more promising area. 
From the theoretical standpoint, this scenario is captured by letting the rv X have a 


multimodal density function, and by relaxing our definition of h(.). Given 1, a positive 





integer, let W" := Uq;,5;),4, <b, <a, <b, <L <a, <b,, and define 


h(u):= ‘max | f fom S06-0)=1} (18) 


a, Sb, Sa, <b, <L <a, <b, ‘ i= 
xeW = 


The constant n represents the number of disjoint sectors that the operator may 


scan before going back to the origin. Clearly, for a, =b,,...,a, =b,, we end up with 
Equation (1) and, moreover, all our developments in this section remain valid with h(.) in 
lieu of h(.). The fact that hi.) is concave (as can be seen by evaluating h"(.) ), implies that 


uh (u)< h(u) for all u = 0, and so strengthens Lemma 1. 


Corollary 5: Suppose 0<7 <0. Then u’ <7. 
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Hil. THE BACK AND FORTH TRAJECTORY 


A. RATE OF DETECTION AS A FUNCTION OF SECTOR LENGTH 

The back and forth trajectory scenario arises for platform mounted sensors, such 
as UAVs and manned aircraft. In this section our goal is to analyze cycle trajectories and 
find an optimal cycle length, where optimal means “maximum expected detection rate 


prior to reneging.” A typical trajectory, with speed 1, is illustrated in Figure 5. 





Figure 5. Sensor location as a function of time for B/F trajectory 


Our setting is the same as for the LTO scenario: X,R, and V(t) are mutually 


independent for all t >0, the speed of the sensor is 1; and the arrival process is Poisson 
with rate 1. With these simplifying assumptions, Figure 5 indicates that because the 
sensor visits the trajectory edges one-half as frequently as the trajectory center, targets 
will renege more often in the edges of the trajectory. Following a similar approach as the 
one used to reach Equation (1), it can be seen that the expected detection rate prior to 


reneging is given by 


sup ‘fr Sy (a+ Hex, u)ds| 


Be) = 2u 





; (19) 


where 
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2(urx) 


c(x,u) = i F r(t)dt + } Fr(tdt,0<x<u. (20) 


The constant = Au) associated with the supremum in Equation (19) is the 
trajectory origin that maximizes the correlation between c(.,u) and the arrival location 
density f,(.). Regarding the function a , It can easily be seen that it is symmetric in 
its first argument around x=u/2 and, taking second partial derivatives with respect to 
x, it is concave non-increasing in x. Figure 6 illustrates Cla. for the exponential 


reneging case. 


c(x) for different values of u for R~exp(10) 


c(x}/(2u) for different values of u for R~exp(10) 
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Figure 6. | Graphic example of C(x,u) with R~Exponential 


The additional reneging present in the B/F trajectory suggests that this policy can 


do no better than LTO; this is made precise in the next lemma. 


Lemma 5. For all u > 0, &u) < g(u). 


Proof outline. The result follows because F,(.)non-increasing implies that 
c(x,u)<2[F, (dt, for all O<x<u. 
We now focus on finding an optimal cycle length = argmax Su). The 


analogue of Lemma 3 to the B/F setting is straightforward. 
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Lemma 6. SupposeP(u<X<i)=1, where w=sup{u:P(X 2u)=1} and 
u = inf {w: P(X <u) =I}, and -0<u<i<oo. Then M<i—u. 


Proof outline. For u, 2u and u_ <u, we have 


- Fy (x)e(x, u)dx [Fe @e(u)de [i Fe @e(au)de 
u,-u) =O) —~*~«S 





(21) 


The other possibilities 


(u_>Uu,u, >U,u, -—u_>u-—U; U_<U,u, <Uu,u,—u_>u—w) are treated similarly. 
Lemma 7 is the equivalent of Lemma 1. 


Lemma 7. Suppose that 0 <7 <0, and that 
(fy(u+ Ao)—7'u,F))u< Fy (u+ &)— Fy (&)—n(u,7), (22) 


for all u>7F.Then M<?7. 


Remark. The last condition is the same as the condition used in Lemma 1, except 


for the error term 


7/2 2x u-r/2 2(u-x) 
} fy(oed| ER } Fy | fused ER } Fo 
n(u,F)=— - : : (23) 





2ER 


The error term becomes insignificant as 7 becomes large, but may be relevant for 


small values of 7 ; see Example 6 below. 
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Proof outline. For u >7 we have 


Fy (ut &) — Fy (&) —9u,7) 





&u) = ER (24) 
u 
and the hypothesis implies that (uw) <Ofor all u>7r. 
Example 6. When P(R=7)=1 for some constant 7 > 0 
c(x,u) = min {(2x) ,7}+min{2(u—x),F}. (25) 


Clearly, &u) = F,(u+ &)—F,,(&%) is non-decreasing for u<7/2. Therefore, if 


the conditions of Lemma 7 apply, 7/2<i%<7 . Once again, the extra reneges that may 


occur at the tails of the B/F trajectory, when compared with LTO, produce %< u' . 


Regarding the general case when both X and R have infinite support, let 


f(u+&) f F,(t)dt 
Pa) = g 





: (26) 
I, f(x+%)c(x,ujdx 4 


Lemma 8. If v'(w)<0 for all w>0, then there exists a unique finite uw’ that 
maximizes g(.) (and is the root of v(u') =0). 


Proof outline: The proof follows along the lines of the proof of Lemma 4, once 


we notice that (uw) = 0 iff Hw) =0.. 
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B. MULTIPLE ARRIVAL CHANNELS 

It is often the case that one UAV is available in a sector of interest with more than 
one concentration of target arrivals. The problem in this case is deciding how many and 
which of these concentrations (referred to here as “channels of arrival”) to patrol in order 
to maximize the overall detection rate. When the channels of arrival are narrow compared 
to the distances between them, they can be approximated by channels of width zero. This 
approximation is valid in various scenarios, such as targets that travel along roads or 
narrow valleys. The main issue in this case is the tradeoff between covering more areas of 
interest and missing detections by “wasting” additional time to travel between the 


channels. 


We will show that this problem can be written as a nonlinear integer 
programming problem. In addition to enabling the solution of the problem by 


enumeration, the formulation provides some insight about these types of problems. 


More specifically, we discuss the case of a sector with n channels of arrival, 


which are all zero in width. We still maintain the assumption that the overall arrival rate 


is 1. The pdf of the arrival locations is defined as P(arrival in channel i) = p, EP = 


i=] 


f. (x) 
P3 
Py DP, 
P2 
----- > 
<r <= oS 
ie ods i, * 


Figure 7. Multiple arrival channels 


Let d,:i=1.n—1 be the distances between the channels, where d, is the distance 


between channel 7 and channel i+1. 
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The patrolled sector length when patrolling from channel k to channel /:/>k 


and back is 


=I 
uy = yd I>k? (27) 


which is also the distance between channels k and /. 


Consequently, the overall detection rate given kK and / is 


a i 
E(k, =— > p c(t, jUz,)- (28) 


Ux) jak 


The problem of deciding which channels to patrol in order to maximize the 


overall detection rate could be written as 


I 
na LY pjclu,,, sf 1<k<l<n. (29) 


tl | QU yy ime 


The inputs to this nonlinear integer programming problem are: 


1. The distances between the channels of arrival 
2. The probability of arrival at each channel 
3. The reneging time distribution (needed to calculate c(u, ;,u,,)). 


Example 7. In the instance of the aforementioned case, targets are expected to 
arrive at one of two roads with identical probability, p =1/2.The distance between the 
roads is d. The dilemma in this case is to choose between surveying one road constantly 
and moving back and forth between the two roads. 


Ze 


In our case, the detection rate can be written as a function of the number of 


channels patrolled m = 1,2 


2) m=1 
Z(m) = } F e(t)dt (30) 


2n-2 
ay 


The above result is very useful in producing a quantitative decision rule. By this 


rule it is beneficial to patrol both roads only if 


[Fe (t)dt 


2 31 
2d GY) 


The interpretation of the above condition is that when the probability of missing a 
target due to traveling between the roads is greater than '4, it is better to survey only one 


road. 


Furthermore, if the reneging time is always shorter that the round trip between the 


two roads, r<2d , this condition becomes ER>d. 


C. CHANGES IN VELOCITY 
We approach the investigation of the effect of the sensor speed in a similar way to 
the LTO case (section 3.3). A change in the sensor’s velocity from 1 to v>Ochanges 


Equation (19) to: 


wolf Fy (xt aje(x/v,u ina 


&,(u) = a (32) 


23 


Changing variables, the latter equals: 


vol (x'+ sea 
a’ \o 


2y 





(33) 


where Sx = Vf 'V) . Because for ¥>9 arbitrary, fx) is a density function, the 
analysis of this section applies to the general velocity setting, with IO in place of 


SQ tt is also evident that “ <“ and g,(u,)<g(u ) if and only if V<!. The meaning 
is that, just as for the LTO trajectory, decreasing the velocity reduces the optimal cycle 


length and the optimal detection rate. 
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IV. INVESTIGATION TIME 


So far we made the assumption that the trajectory of the sensor is predetermined 
and does not change due to real time events. In certain situations, however, targets are 
investigated by the sensor upon detection, during which time the sensor remains 
immobile. This type of sensor policy is common when the identification of the target is 
required before a decision is made as to the necessary reaction, or when the sensor has to 
“escort” the target until the reaction force arrives at the scene. In these cases, a dynamic 


real time algorithm is needed in order to maximize the MOE. 


Specifically, suppose that each target is investigated « units of time after 
detection, and that 7(x,t) represents the clock when position x was last visited. In the 
LTO setting, we assume that at any instant the camera operator has the option of 
continuing to move rightward, or it can make an instantaneous leap to any leftward point. 
Our MOE is “expected time until the first detection,” and we wish to select a leftward 
leaping point or continue moving rightward in order to minimize that MOE. In the B/F 
framework there are only two options: move towards the left or to the right, so the 
analysis is somewhat simpler. In either case, we make the same basic assumptions: the 
rv’s X and R are independent, and independent of the arrival process, which is 
homogeneous Poisson; we also assume that sensor velocity is 1. Clearly the setup is 
markovian, in that future detections depend only on the current sensor position and on the 


latest visiting times to each location in the perimeter. Let z(¢,) be the sensor’s position at 


time, , and suppose the sensor does not change course (leap in LTO or change course in 


B/F). Then 


. : . 7 ot x- (H(t +x) = 
P(no detections in next At units of time) = I ‘i Ft (a(t)) +x) exp(— I Ne Oe ? (2)dz Ar ; 


which follows because for each position x, the detection process is Poisson with rate 
F,(t),0 <t < (x) (cf. the discussion preceding Equation (2)), and velocity is 1. The last 
expression can be used to compute the expected time until the first detection without 


changing course, D (7,7), given position 7, last visiting clock 7 


no change 
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Dayo change (4577) = [: P(first detection occurs after t)dt 


n 


= iN P(no detection until t)d¢ . (34) 


= A [E JS (a(t)) +x) exp (-(re F, (2)dz dd 


0 
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V. FURTHER STUDY 


The cases studied in this paper assume that the target arrival process is a 
homogenous Poisson process, and that reneging times and arrival locations are 
independent of each other and of the Poisson process. An important case is when the 
reneging times are dependent on the targets’ arrival locations. Such target behavior can 
occur when varying terrain conditions in different locations along the patrolled perimeter 
affect the crossing times of the targets. As explained in the Introduction, the detection 
process remains Poisson under conditional independence, so that the approach taken in 
this paper to find optimal cycle lengths remains valid. Further work is required to explore 


a dependence structure that keeps our formulation analytically tractable. 


An important class of scenarios arises when multiple sensors can be used in a 
sector. Analysis of this case boils down to the derivation of optimal employment of more 
than one sensor, the system’s performance, the cost structure, and sensitivity to the 
number of sensors. Equivalently, there may exist multiple target types, with a reward 


associated with each type; the goal here would be to maximize the reward earning rate. 


Regarding the target arrival process, arrivals in different time periods may be 
correlated, which destroys the Poisson assumption. A more realistic model would allow 
for these dependencies, also possibly taking into account spatial effects. Markov random 


fields, whose theory is well developed, seems the ideal candidate to explore this venue. 


In our analysis it is assumed that the sensor’s location does not affect the targets 
behavior; this assumption can be made for hidden sensors like hidden cameras and high 
altitude UAV’s, but might not always hold. This suggests a game-theoretic facet to our 
problem. Specifically, targets may be smart and behave in order to lessen the probability 


of capture. 
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VI. CONCLUSIONS 


This paper presents a model for a moving sensor patrolling a perimeter. The goal 
is to find a patrol policy for the sensor that minimizes the infiltration across the perimeter 
(maximize detections). We assume the sensor follows a periodic trajectory at constant 
speed. The controlled parameters (decision variables) of the trajectory are the beginning 


and end points of the patrolled sector. 


Theoretical framework is established for two general trajectories: a Back and 


forth (B/F) trajectory and a Leap to origin (LTO) trajectory. 


The LTO trajectory represents cases where the perimeter is scanned in a constant 
speed (or constant angular speed in the case of a stationary camera) from the origin to a 


certain point, and then instantaneously starts to be scanned from the origin again. 


The target behavior is described extensively in the introduction. The two main 
assumptions in this paper are that the target arrival process is a constant-rate Poisson 


process and the target arrival locations, arrival times and reneging time are independent. 


The assumption that the sensors speed is | is made for convenience. Section 3.3 


discusses the adaptation of the results for general speed. 


For the LTO trajectory, we derived an expression for the rate of detection given 
the reneging, arrival location distribution and the length of the patrolled sector (the 
decision variable). In Lemma | it is shown that it is never beneficial to patrol a sector 
larger than the upper bound on the reneging time distribution, if outside it the pdf of 
target arrival location is non-increasing. 

In a case where the variability in reneging times is very low, the reneging time 


can be represented as deterministic and r is the optimal sector length. 


An interesting private case is when the arrival location distribution is uniform and 


r is smaller than the width of the arrival distribution, u. In this case, the loss in detection 
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rate due to the increase in the sensor revisit times is equal to the gain in detection rate due 
to the increase in covered sector for all u €[r,u]. The result is that all the sector lengths 


in that interval are optimal. 


We then consider a different scenario. In this scenario, targets appear on a 
perimeter and do not renege. This scenario is relevant to cases where a crossing 
malevolent agent leaves a track at the location of its infiltration (for example, tracks on a 
trail). The objective in these situations is to detect the tracks left by the target within the 
shortest possible time. The sensor platform in this case could be a motorized ground 
patrol or a visual scanning sensor. Example 2.b shows that maximizing detection rate in 
the case of uniform reneging time is equivalent to minimizing the expected time to 


detection of targets (waiting time) in the aforementioned scenario. 


A basic limitation of the suggested model is that the sensor footprint is modeled 
as a one-dimensional point. This assumption is a reasonable representation of reality 
when the size of the sector is significantly greater than the sensor’s footprint width. 
Another interpretation of the model is that it represents only the component of detection 
that is produced by the movement of the sensor. For example, a stationary sensor has a 


detection rate of zero in our model, while in reality this is not necessarily true. 


For the B/F trajectory, we also derived an expression for the rate of detection 
given the reneging, arrival location distribution and the length of the patrolled sector. In 
this case, unlike the LTO trajectory, the probability of target detection given its arrival 
location along the patrolled sector is not necessarily the same for all locations. In fact it is 
not the same for most realistic scenarios. This fact is captured by the expression c(x,u), 
which represents that conditioned probability: P(detection|u,x:x €[0,u]) =c(x,u). 


This function is determined by the sector length and the reneging distribution. 


A very close relationship exists between the detection rate in the B/F trajectory 


and the LTO trajectory and it is analyzed in the B/F section (IID) of this paper. 
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